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A detailed numerical and analytical examination of the evolution of stochastic magnetic fields 
between a putative magnetogenesis era at high cosmic temperatures T ~ 100 MeV - 100 GeV and the 
present epoch is presented. The analysis includes all relevant dissipation processes, such as neutrino- 
and photon- induced fluid viscosities as well as ambipolar- and hydrogen- diffusion. A simple and 
intuitive analytical model matching the results of the three-dimensional MHD simulations allows for 
the prediction of pre-recombination and present day magnetic field correlation lengths and energy 
densities as a function of initial magnetic field energy density, helicity, and spectral index. Our 
conclusions are multi fold, (a) Initial primordial fields with only a small amount of helicity are 
evolving into maximally helical fields at the present. Furthermore, the simulations show a self- 
similarity in the evolution of maximally helical fields implying a seemingly acausual amplification of 
magnetic fields on large scales is observed, (b) There exists a correlation between the strength of the 
magnetic field B at the peak of it's spectrum and the location of the peak, given at the present epoch 
by: B w 5 x 10 -12 Gauss (L/kpc), where L is the magnetic field correlation length determined by 
the initial properties of the magnetic field, (c) Concerning studies of generation of cosmic microwave 
background (CMBR) anisotropics due to primordial mag netic fields of B ~ 10" 9 Gauss on > 10 Mpc 
scales, such fields are not only impossible to generate in early causal magnetogenesis scenarios but 
also seemingly ruled out by distortions of the CMBR spectrum due to magnetic field dissipation on 
smaller scales and the overproduction of cluster magnetic fields, (d) The most promising detection 
possibility of CMBR distortions due to primordial magnetic fields may be on much smaller scales at 
higher multipoles I ~ 10 6 where the signal is predicted to be the strongest (e) It seems possible that 
magnetic fields in clusters of galaxies are entirely of primordial origin, without invoking dynamo 
amplification. Such fields would be of (pre-collapse) strength 10~ 12 — 10~ n Gauss with correlation 
lengths in the kpc range, and would also exist in voids of galaxies. 

PACS numbers: 

I. INTRODUCTION 

Magnetic fields exist throughout the observable Universe. They exist in stars, in the interstellar medium, in galaxies, 
and clusters of galaxies (for a review see 0]), where in the latter two environments they are often observed with /xGauss 
strength. Magnetic fields likely also reside in the intergalactic medium, though at present, their strength may only be 
limited by Faraday rotation measures of distant quasars 2] . The origin of galactic- and cluster- magnetic fields is still 
unknown. A plausible, though by far not convincingly established, possibility is the generation of magnetic seed fields 
and their subsequent amplification via a galactic dynamo mechanism. Seed fields may be due to a variety of processes 
(and with a variety of strengths), such as the Biermann battery within intergalactic shocks Q, stellar magnetic fields 
expelled in planetary nebulae, or during supernovae explosions, either into the intragalactic, or in the presence of 
galactic outflows into the intergalactic medium , as well as due to quasar outflows of magnetized plasma Q . Seed 
fields may also be of primordial origin with a multitude of proposed scenarios. These include generation during 
first-order phase transitions (e.g. QCD or electroweak) , around cosmic defects, or during an inflationary epoch (with, 
nevertheless, extremely small amplitudes), as well as before the epoch of neutrino decoupling or recombination. For 
a review of proposed scenarios we refer the reader to [y, |jj . 

The philosophy in prior studies of primordial magnetogenesis is often (but not always) as follows. After establishing 
a battery mechanism (e.g. separation of charges and production of currents) and a "prescription" or estimate for the 
final, non- linearly evolved magnetic field strength (e.g. equipartition with turbulent flows), subsequent evolution is 
approximated by simply assuming frozen-in magnetic field lines into the plasma. Though this may be appropriate on 
the very largest scales, it should be clear, that this may not be the case on the fundamental coherence scale of the 
field. Here, coupling of the magnetic fields to the gas induces non-linear cascades of energy in Fourier space. The 
characteristics of initially created magnetic field are thus vastly modified during cosmic evolution between the epoch 
of magnetogenesis and the present. The final step in such studies is then often to determine field strengths on some 
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prescribed scale (e.g. 10 Mpc) typically falling in the range 10 -30 Gauss < B < 10~ 20 Gauss, inferring that this may 
be sufficient to seed an efficient dynamo for the production of galactic- and cluster- magnetic fields. This is observed 
in negligence of the fact that much stronger fields on smaller scales, result not only from a variety of astrophysical 
seeds, but from these very same primordial scenarios. 

Considering the likelihood of a magnetized early Universe (i.e. due to the large number of charged particles 
and the multitude possible of out- of- equilibrium processes) it should be instructive to develop a somewhat complete 
picture of magnetic field evolution in the early Universe, subsequent to the epoch of magnctogcncsis. This should be 
accomplished irrespective of such fields providing the seeds for galactic fields, or not. For example, it may be that 
at some later time relatively weak field strength in galactic voids are measurable via the propagation of the highest 
energy- cosmic rays 0, E| , or via accurate measurements of 7-ray bursts 0] . The interpretation of such putative 
measurements, which could hint to fields of primordial origin, is then possible only if one understands the evolution 
of these fields between magnetogenesis and the present. 

One step in this direction has been performed by Dimopoulos & Davis 11] as well as Son [l2j who exactly considered 
such non-linear processing of magnetic fields due to magnetohydrodynamic cascades in the early Universe. Another 
step has been provided by Jedamzik et al. 0] (hereafter; JK098) and shortly after Subramanian & Barrow 0] 
who considered fluid- viscosity (due to neutrinos and photons) induced damping of magnetic fields. The study by 
Son, though describing appropriately the gross non-linear features of MHD evolution, does not properly deal with 
the effect of fluid viscosity [53j, with the net effect of estimates of present day coherence lengths being orders of 
magnitude smaller than those we find. Moreover, this study, as most others, does not provide explicit expressions for 
the final magnetic field energy. The study by Ref. on the other hand, though examining the effect of photon drag 
before recombination, employs a somewhat particular model of magnetic field coherence length growth which is not 
supported by results of numerical simulations. The study by JK098 is strictly speaking applicable only in the linear 
regime (i.e. under the assumption of a homogeneous background magnetic field), whereas Subramanian & Barrow also 
considered a limited class of non-linear configurations. Results of both studies on the magnetic field coherence length 
at recombination are identical to those found in the non-linear analysis attempted here. Nevertheless, at intermediate 
stages of evolution (i.e. well before recombination) the predicted magnetic coherence length in these studies deviates 
from that found here. Moreover, none of these works discusses the effects of ambipolar diffusion after the epoch of 
recombination, neither verifies claims by complete three-dimensional numerical simulations. 

Two-dimensional- (e.g. |l5)) and three-dimensional- (e.g. |l6j) numerical simulations of magnetohydrodynamics in 
the early Universe were performed in the context of maximal helical fields. Similarly, effective three-dimensional 
cascade models [l2T ] have also been emplo yed . Helicity of primordial magnetic fields could play an important role, 
as noted by a number of authors 0, 0> lisl 0, H3, , as it may significantly speed up the growth of magnetic 
field coherence length, thereby leading to potentially large magnetic fields strengths on comparatively large scales 
(~ 1 — lOkpc depending on the amount of initial helicity). It has also been argued that a net primordial magnetic 
helicity may be potentially linked to the cosmic baryon-to-entropy ratio (e.g. 0, H3|). Adopted models of field 
evolution are either appropriate to turbulent evolution |l9j . or to viscous evolution (i.e. assuming a drag force due to 
photons [HI)- Before passing, we also note studies of the effects of magnetic fields on the cosmic microwave background 
radiation (CMBR), as for example, the generation of temperature anisotropies below [22j,|23j or above 1241 the Silk 
damping scale, as well as the distortions of the CMBR Planck spectrum by magnetic field dissipation [25j. These 
studies have also to include certain evolutionary features of magnetic fields, but due to the largeness of the Silk 
scale (~ 10 Mpc comoving) back reaction of the peculiar flows generated by magnetic fields on the magnetic fields 
themselves. This is in contrast to the im por tance of back reaction on the typically much smaller magnetic field 
coherence scale (such as for the analysis in |25|) . 

In this paper we attempt to provide a unified picture for the gross features of magnetic field evolution in the early 
Universe. As a function of initial conditions for the magnetic fields generated during a putative magnetogenesis era, 
we predict the magnetic field coherence length and magnetic energy density for all subsequent epochs for fields of 
arbitrary strength and helicity. Our treatment incorporates all the relevant dissipative processes, in particular, due 
to photon- and neutrino- diffusion as well as free-streaming, and due to ambipolar- and hydrogen- diffusion. 

The outline of the paper is as follows. While many of the preliminaries to the discussion, such as the equations, 
treatment of Hubble expansion, and magnitudes of dissipative terms, are deferred to the appendices, Sec. [n] immedi- 
ately commences with a discussion of turbulent MHD cascades and the presentation of results of three-dimensional 
numerical simulations. In Sec. IIIII magnetic field evolution in the viscous regime before recombination (which is a 
regime particular to the early Universe) is discussed and numerically simulated, whereas Sec. IIVI discusses the effects 
of ambipolar diffusion after recombination. The general picture and detailed analytical results for cosmic magnetic 
field evolution is developed in Sec. whereas Sec. IVII provides a discussion of the highlights of our findings. In 
the Appendices 1X1 and iBl we compile the MHD equations appropriate for the study of magnetohydrodynamics in the 
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expanding Universe, whereas Appendix IO complies the various dissipation terms in the early Universe. Details on the 
generation of helical fields are given in Appendix [D] and details on the numerical simulations in Appendix lEl 

II. TURBULENT MAGNETOHYDRODYNAMICS 

In this section we discuss general features of the evolution of magnetized fluids in the turbulent regime (Reynolds 
number R e 3> 1 as applicable well before neutrino coupling and recombination) , such as the decay of energy density 
as well as the growth of magnetic field coherence length. The exceedingly large Prandtl numbers (cf. Appendix |0 
in the early Universe allow one to neglect dissipative effects due to finite conductivity. Further, the generation of 
primordial magnetic fields in magnetogenesis scenarios is generally believed to occur during well-defined periods (e.g. 
QCD-transition). Subsequent evolution of these magnetic fields is therefore described as a free decay without any 
further input of kinetic or magnetic energy, i.e. as freely decaying MHD. Due to the largeness of the speed of sound in 
a relativistic plasma v s = I/v 7 ^, the assumption of incompressibility of the fluid is appropriate during most epochs, as 
well as for a large range of initial magnetic field configurations and energy densities. Exception to the incompressibility 
may occur for initial conditions which result in magnetic fields of strength B > 6 x 10 -11 Gauss (comoving to the 
present epoch, cf. Sec. and only after the decoupling of photons from the flow. 

To verify theoretical expectations we have performed numerical simulations of incompressible, freely decaying, 
ideal, but viscous MHD. These simulations are performed with the help of a modified version of the code ZEUS- 
3D |26l 1271 128 | in a non-expanding (Minkowski) background. Modifications lie in the inclusion of fluid viscosities, e.g 
a drag coefficient a as given in Eqs. ©, 1C8(1 . and l|C9(l . From the discussion in Appendix 151 it should be clear, that 
for most purposes results of numerical simulations with existing (or slightly extended) codes with Minkowski metric, 
may be re-interpreted into results of MHD in an expanding Universe with FRW metric, when rescaled variables as 
given in the Appendix are considered. For details of the numerical simulations the reader is referred to Appendix [E] 

Incompressible MHD is described by the following equations (for the equations of compressible MHD the reader is 
referred to Appendix lAl and ITj|) : 



dt 

dv A 



" V +(vV)v-(v A -V)v A = f , (1) 



dt 



+ (v V)v A - (v A - V)v = iA7 J v A , (2) 



where we have defined a local Alfven velocity va(x) = B(x)/ y/ Att{q + p), and where v, B, g and p are the velocity, 
magnetic field, mass-energy density, and pressure, respectively. Here fluid dissipative terms in the Euler equation are 
given by 

f ?yV 2 v l mtp < I 

f=\ > (3) 

( -av l m f P > I 

where there exists a distinction between dissipation due to diffusing particles, with mean free path smaller than the 
characteristic scale Z m f p -C L, or dissipation due to a free-streaming (i.e. / m f p 3> L) background component exerting 
drag on the fluid by occasional scatterings with fluid particles. Both regimes are of importance in the early Universe 
as already noted in JK098. Note that in the computation of va only those particles with Z m f p <C L contribute to g 
and p. An important characteristics of the fluid flow is given by it's local kinetic Reynolds number 

^rnfp "C I 

mn = ~ = { 11 , (4) 

Indp 3> I 

with I some length scale. The Reynolds number is a measure of the relative importance of fluid advective terms and 
dissipative terms in the Euler equation, given by the ratio of a typical dissipative time scale — (I 2 /r], 1/a) to the eddy 
turnover time scale r e ddy — l/v. For most magnetic field configurations it is possible to define an integral scale, L, i.e. 
the scale which contains most of the magnetic- and fluid kinetic- energy. We will frequently refer to this scale as the 
coherence scale or coherence length of the magnetic field. In the case of turbulent flow, with R e (L) I on this scale, 
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the decay rate of the total energy is independent of dissipative terms and only depends on the flow properties on the 
integral scale. This is in contrast to the decay of magnetic- and fluid- energy in the viscous regime, R e (L) <C 1, where 
the total decay rate depends on the magnitude of viscosities. In the following, the dynamic evolution of magnetic 
fields in the former (turbulent) regime will be studied. 



A. Nonhelical fields 



Consider Eqs. JIJ and J5J) with a stochastic, statistically isotropic, magnetic field and, for the purpose of illustration, 
with initially zero fluid velocities. For the moment we will also assume that the magnetic field does not possess any net 
helicity. In the limit of large Reynolds numbers on the coherence scale, the dissipative term may be neglected on this 
scale. Magnetic stresses (the third term on the LHS of Eq. JJJ) will establish fluid motions of the order v ss va within 
an Alfven crossing time ta ~ 1/va, at which point back reaction of the fluid flow on the magnetic fields will prevent 
further conversion of magnetic field energy into kinetic energy. The resultant fully turbulent state is characterized by 
close-to-perfect equipartition (in the absence of net helicity) 



,2\ ^ /., 2 



v- 



va 2 ) , (5) 



between magnetic and kinetic energy. This may be seen in Fig. ^ which shows the decay of magnetic- and kinetic- 
energy in freely decaying turbulent MHD. 

Non- linear MHD processes quickly establish turbulence on scales below the integral scale (cf. Fig. Working 
with Fourier transforms (assuming statistical isotropy and homogeneity) and defining the total magnetic- and kinetic- 
energy density 



E w / d In k k s 



((M 2 > + (Kfcl 2 )) = f dlnfc^, (6) 



one finds that a typical root-mean-square velocity perturbation on scale I — 2n/k is vi y/k 3 (\vk\ 2 ) ~ \[Ei. Note 
that in the above and for the remainder of this section we set (g+p)/2 = 1 (cf. Eq. £3), as frequently done in studies 
of incompressible MHD, such that energy density has the dimension of velocity square. By inspection of the Fourier 
transformed Eqs. I|C3(1 and l|C4(l it may be seen that dissipation of energy is dominated by flows on the smallest scales 
(largest k), given that energy spectra E\ fall not too steeply with growing k. Dissipation of energy into heat thus 
occurs at some much smaller scale Z^iss <C L (where R e (l) « I). The transport of the fluid energy from the integral 
scale L to the dissipation scale Zdi ss occurs via a cascading of energy from large scales to small scales, referred to as 
direct cascade. 

Ever since the work of Kolmogorov, as well as Iroshnikov and Kraichnan it is known that this cascading of energy 
occurs as a quasi-local process in fc-space, with flow eddies on a particular scale I breaking up into eddies of somewhat 
smaller scale ^1/2. This continuous flow of energy through k-space 

« — ~ const(fc) , (7) 
at ti 

results in a quasi-stationary energy spectrum on scales I < L, with energy flow rates approximately independent of 
wave vector. We remind the reader that throughout L denotes the integral scale. Typical dissipation times 77 are 
given by an eddy-turnover time scale t\ ~ T e ddy,z ~ l/vi ~ l/(>/E~i) in the unmagnetized case, and 77 ~ tjk,1 ~ 
(1/vi)(va.l/ v i) ^ I /{Ei) in the MHD case. Here the latter time scale takes account of the Alfven effect, in particular, 
a prolonging of energy cascading time in the case of scattering of oppositely traveling Alfven waves [2^. l30j |. A 
description of the turbulence as a superposition of Alfven waves seems appropriate, as the magnetic field on the 
integral scale («a z, 3> v A,i) effectively provides a strong and locally homogeneous background field on the scale 
I <C L. Note that, T c( jdy,L = 7"iK,l on the integral scale itself. 
Equation Q may thus be solved yielding 

/c~ 5//3 : unmagnetized 

/•:, /•-■/ /.■ x { , (8) 

k~ z / 2 : magnetized 



where Ej- = Ei, to yield the familiar Kolmogorov- and Iroshnikov-Kraichnan (IK) spectra on scales below the integral 
scale. It is matter of current debate if spectra as proposed by IK indeed result during MHD turbulence. Goldreich & 
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Sridhar |3l| have established that MHD turbulence is intrinsically anisotropic with eddies elongated in the direction 
of the background (integral scale) magnetic field (i.e. ku ~ fc^ 3 L -1 / 3 where ku and k± are wave vectors parallel 
and perpendicular to the background magnetic field «a l) an d energy cascading more rapidly in fcj_-space orthogonal 
to the magnetic field. Though the predicted anisotropy has been observed in numerical simulations, the predicted 

~ —5/3 

modification of IK spectra, in particular, the existence of "one-dimensional" Kolmogorov-type spectra Ek oc k , 
has been not |32j. Rather, these spectra seem to follow the one proposed by IK. In contrast, Miiller & Biskamp^33| 
find "three-dimensional" energy spectra consistent with Kolmogorov but inconsistent with IK. Inspection of Fig. [21 
which shows such "three-dimensional" spectra in our 256 3 numerical simulations of freely decaying MHD, illustrates 
how difficult it is to distinguish between exponents —5/3 and —3/2. This is due to the inertial range between the 
dissipation scale (here given by a few times the Nyquist frequency due to numerical dissipation) and the integral scale 
being rather small. Moreover, both scales do not seem to be well defined, resulting in a small-scale spectrum more 
consistent with an exponential than a power-law. This is also not significantly changed when one proceeds to 5I2 3 
simulations, such that a numerical confirmation of one, or the other, spectrum may be premature. 

Evolution of global properties of the magnetic field in freely decaying MHD, such as total energy density and 
coherence length, depend on the magnetic field spectra on scales above the integral scale, I > L, and are related to 
initial conditions. Consider an initial magnetic field with blue spectrum, 

Ek ~ E ° = E ° (i) for 1 > L ° ■ (9) 

The scale-dependent relaxation time, t; ~ 1/va.i ~ 1/yEi (with va.i = \/k 3 (\vA.k\ 2 }) increases with scale as t; oc 
l 1+n / 2 . Transfer of magnetic energy to kinetic energy and a fully developed turbulent state may only occur for times 
t > ti . When such a state is reached the energy on this scale decays through the cascading of large-scale eddies to 
smaller-scale eddies down to the dissipation scale. Since the relaxation time for the "next" larger scale I is longer, 
this larger scale now becomes to dominate the energy density, i.e. becomes the integral- or coherence- scale. This is 
sometimes referred to as selective decay of modes in k-space. The remaining energy density is then the initial energy 
density of modes between the very largest scales and this next larger scale. Given these arguments and the initial 
spectrum of Eq. @ one then may derive for the time evolution of energy and coherence length of the magnetic field 

EkE o(J-) L ~ L o {yY^ no helicity, R e > 1 , (10) 

for t > To, where To is the relaxation time on the scale Lq, i.e. To ~ Lq/^/Eq ~ Lq/v^ , and where indices denote 
quantities at the initial time. For instance, for a spectral index of n = 3 (which corresponds to the large-scale magnetic 
field due to a large number of randomly oriented and homogeneously distributed magnetic dipoles |34j) the energy 
density follows E oc i~ 6 / 5 which is Saffman's law known from fluid dynamics |35l l36| . 

An increase of magnetic field coherence scale with time due to selective decay may be observed in Fig. [21 whereas 
the decay of magnetic energy density for a variety of initial magnetic field spectra is shown in Fig. [3] It can be 
seen that initial spectra with larger n indeed lead to a more rapid decrease of energy with time as predicted by 
Eq. 1|10|) . Nevertheless, comparison of the theoretically expected decay exponents (cf. Eq. (|f 0\i ) to the numerically 
found exponents (cf. Fig. [3J indicate slight differences. Generally, our numerical simulations result in a slower energy 
decay than predicted by Eq. (|1U|) . For example, the theoretically predicted damping exponent, i.e. E(t) oc i -7 , for 
a n = 3 initial energy spectrum is 7 = 1.2 whereas the best fit of our numerical simulations gives 7 « 1.05. It is 
not easy to find a physical explanation for this, as it would entail an additional with time increasing slow-down of 
relaxation at large scales I > L and/or slow-down of energy dissipation of already turbulent modes at small scales 
I < L. In either case, to explain such a phenomenon a quantity with physical dimension of length or velocity, which 
has not yet entered the analysis, should exist. Given that the assumed initial magnetic field distribution is statistically 
self-similar on different scales, and that helicity is negligible, this quantity may only be the dissipation length and/or 
length of the simulation box. Whereas the latter is a complete numerical artifact, the former is so widely separated 
from the integral scale during most periods of the high Reynolds number flow in the early Universe, that we expect it 
not to influence the dynamics on the integral scale. We have noted, that spectra at late times show a peak region AL 
quite spread, and are likely only marginally resolved by the simulations. In any case, larger numerical simulations are 
required to address this effect seen also by others (e.g. [lH I37I 1. 
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FIG. 1: Comparison of the time evolution of the magnetic (solid line) and the kinetic (dashed line) energy in the turbulent 
regime (R e 3> 1) for a magnetic field without initial helicity. For comparison, also the theoretical damping law, E oc t~ 1,3 , is 
shown (dotted line). Here, the simulation was performed on a mesh with 128 3 grid points and the magnetic field is excited up 
to k c « f6 with a spectral index n ~ 4 (cf. Eqs. JUJ and (llUl and Appendix El. 

B. Helical fields 

We have so far studied the evolution of a statistically isotropic and homogeneous magnetic field in the absence of 
net helicity (see appendices IU1 and iDl for the definition and dissipation of magnetic helicity). Given that magnetic 
helicity should be an ideal invariant in the early Universe, and that magnetic fields with even small initial net helicity 
ultimately reach maximal helicity density 

U < H max w (B 2 L) « (8tt) EL, (11) 

it should be of importance to also study the maximally helical case. Note that, a maximally helical state is reached 
during the course of MHD turbulent evolution due to a slower decay of the helical component of fields as compared 
to the non- helical one (cf. Eq. i|l(J|) and Eq. (|13f) below). When maximal helicity is reached magnetic field evolution 
is significantly altered with respect to the case of zero, or sub-maximal helicity. Fig. 0] shows the results of 256 3 
simulations of the evolution of the ratio between kinetic- and magnetic- energy density Y assuming initial conditions 
of a maximally helical field and negligible velocity perturbations. After a relaxation time of the order of the Alfven 
crossing time L/v\ L over the integral scale a quasi-steady state with constant T w 0.2 develops. Note that in contrast 
to the turbulent, non-helical case, full cquipartition is not reached. The associated spectrum of T is shown in Fig. [5] 
showing that at the integral scale kinetic energy density is always smaller than magnetic energy density. Though 
not apparent from the figure, an integral scale (i.e. the scale of maximum energy density) defined for kinetic energy 
density only L km , trails the integral scale for total energy density with time, i.e. L km /L < 1 with a ratio approximately 
constant in time. Magnetic field spectra for this simulation are shown in Fig. EI It is seen that inertia! range magnetic 
spectra at I < L are well-described by power-laws over a limited range in k-space. For this exponent (n = 4) we find 
Ek oc fc' 3 , [3 w —1.7, significantly steeper than either Komogorov or IK. 

Fig. also illustrates the intriguing property of self-similarity of spectra at different times. This phenomenon of 
self-similarity has also been observed by [16j. Magnetic field amplification on very large scales occurs even at times 
much shorter than the typical relaxation time for magnetic fields (i.e. Alfven crossing time) on these scales, indicating 
the topological constraint (by helicity) imposed on the field evolution. Note that if magnetic fields on large scales 
would not be enhanced, magnetic coherence length could not grow with time, as generally the initially existing energy 
density on large scales would not suffice to keep H constant. Having performed simulations of maximally helical fields 
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FIG. 2: Evolution of magnetic energy spectra in the turbulent regime for a magnetic field with no initial helicity. Here, the 
spectral index of the initial energy spectra is n w 4. Note that Ek as opposed to Ek is shown (cf. Eq. JSJ 



with different initial spectral indices n we have noted that though the amplitude of large-scale magnetic field grows 
with time, the spectral index of the magnetic field configuration on large scales seems to be approximately preserved. 

There seems to be a misconception in the literature (see, e.g. [l9Ll20j ') that maximally helical fields do not dissipate 
energy via excitation of fluid flows and the subsequent dissipation of these flows due to fluid viscosities. It is argued, 
that maximally helical fields with a fairly peaked spectrum are essentially force-free (i.e. va x (V x va) ~ 0) and 
may thus not excite fluid flows. Note that if this indeed would be the case, Eq. would imply trivial magnetic field 
evolution B = const for initially zero velocity fluctuations and resistivity. Though the magnetic stresses in the Euler 
equation are indeed smaller for a maximally helical field as compared to a non-helical field of similar strength, an 
increase of magnetic coherence length and the continuous excitation of sub-equipartition fluid flows are observed in 
our simulations. In the limit of large Prandtl number, dissipation of these flows by fluid viscosity will then provide 
the main dissipation of magnetic field energy. 

The decay rate of total magnetic energy in freely decaying MHD turbulence of maximally helical fields may be well 
approximated by the decay rate of energy on the integral scale 



dE 
~dt 



E 

TL 



£3/2 



-r 



£5/2 



-v. 



(12) 



where tl ~ L/va l, and Eq. (|llfl for a maximally helical field has been employed in the second step. Since Ti, and 
r (see Appendix El and Fig.^J) are constant it is straight forward to derive the power-law exponents for the decay of 
energy and growth of coherence length with time 



E ~ En ( — 

To 



-2/3 



Lft ( - 



2/3 



maximal helicity , R e 3> 1 , 



(13) 



for t > To ~ Lq I ' \fEft~ ss Lft/v^ , yielding a predicted decay which is independent of the spectral index of the large-scale 
magnetic field. The correctness of Eq. (|13f) has been recently questioned by Biskamp & Miiller [38j]. These authors 
advocate a decay of kinetic energy with time as Y cx E/H, yielding a modified Eq. (O dE/dt ~ E 3 /H 3/2 , and energy 
decay E cx t -1 / 2 . We note here that a decay of T was not found in our simulations. Moreover, a relationship T cx E/H. 
is dimensionally incorrect, and must be modified by an as yet unknown quantity of dimension length. Due to the 
absence of a physically well- motivated choice for this quantity (other than Z^iss or £box), we suspect their results to be 
an artifact of limited resolution. In particular, Biskamp & Miiller |38| observe a decay in T only at late times, when 
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FIG. 3: The evolution of the magnetic energy in the turbulent regime for different initial energy spectra n, where Ek = 
k 3 |fofc| 2 oc k n with a cut-off k c « 32. Here, the initial magnetic field is non-helical. In this case, the damping law depends on 
the spectral index (cf. Eq. 1101 ). For comparison, the theoretical predicted damping laws for n = 1 (E oc t~ 6r ) and for n — 5 
(E oc i -1 ' 4 ) are also shown. 



the coherence scale has already moved dangerously close to Lbox- Note, that larger kinetic (numerical) viscosities 
result in larger magnetic dissipation times (cf. Sec. 1111(1 . Therefore, the rather moderate Reynolds numbers (0(1O 3 )) 
achieved in numerical simulations could be responsible for the slower decay rates found in these simulations. 

Fig. [7| shows the total magnetic energy as a function of time for a variety of maximally helical magnetic fields of 
different initial spectral index. With the exception of the rather red spectrum n = 1, for which the Fourier transform 
of helicity is not peaked in k-space, the decay of energy seems to be indeed approximately independent of spectral 
index. Residual dependencies on n may possibly be associated with the non-conversation of helicity as shown in 
Fig El This dissipation of helicity in our simulations is due to numerical dissipation at the Nyquist frequency. Similar 
to the case of non-helical fields, the decay slopes observed in the simulations are somewhat shallower than those 
predicted by Eq. I|13[> . For example, for n = 5 we find a damping exponent of 7 w 0.5 (coincidentally agreeing with 
|38|). Arguments very similar to those presented at the end of the previous paragraph, in particular, the absence of 
a quantity of dimension length or velocity beyond those employed in Eq. l(T3)l , makes us believe this deviation to be 
unphysical. 



III. VISCOUS MAGNETOHYDRODYNAMICS 



Magnetic field dissipation in high Prandtl number fluids may also occur in the viscous regime, where kinetic 
Reynolds numbers are much smaller than unity. Of particular importance to MHD evolution in the early Universe is 
the case of photons or neutrinos free-streaming over the scales of interest, l m f p 3> I, resulting in a drag force in Eq. 
with drag coefficient a. To our knowledge, such a case of "dragged" MHD is not encountered in other astrophysical 
environments, and numerical simulations of this case have so far not been performed. 

Consider again the Eulcr equation Eq. Whereas in the turbulent case there is a balance of the terms on the 
LHS, which are all of similar magnitude, in the dragged case there is a balance between the last term on the right- 
hand-side and the dissipative term /, with all other terms negligible. In the terminal velocity regime one finds thus 
(using Eq. ©) 

v « - (v A • V) v A (14) 
a 
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FIG. 4: Time evolution of T = -Bkin/^mag for maximal helical magnetic fields with different spectral indices n in the turbulent 



and magnetic energy is not established for helical magnetic fields. 



such that vi w i>a,l (Tdra.g/TA.,0 "C va.i for Tdrag = ol 1 <C ta,i- This yields a kinetic Reynolds number of 



*e«(^)«l. (15) 

Though one would naively expect that at small Reynolds number the total energy gets immediately dissipated due 
to viscous terms, this is not the case (JK098). For large Prandtl number the energy may only be dissipated via 
the excitation of fluid motions. Nevertheless, due to the strong drag, such excitation is slow and inefficient, and a 
system with r<C 1, i.e. well below equipartition between magnetic- and kinetic- energy results. Since the dissipation 
rate is proportional to the velocity fluctuations v the net effect of strong fluid viscosities is a delayed dissipation and 
quasi-frozen-in magnetic fields. Note that in the case of viscous MHD, flows are effectively dissipated on the integral 
scale, and cascading of energy in fc-space is not required. One finds for the energy dissipation rate 

"757 ~ TT~ > Wlth TL ~ L Vl W ' ( 16 ) 

at tl L z a E 

and where tl a formal eddy turn-over time scale identical to the overdamped time scale for the evolution of overdamped 
Alfven and slow- magnetosonic modes as found in JK098. 



2 



A. Nonhelical fields 



With a blue spectrum (n > 0) for the magnetic fields on large scales as given in Eq. @, and with very similar 
reasoning as in the turbulent case, one may compute the asymptotic power-law for decay of energy density and growth 
of magnetic field coherence length as 



E « En 



t 



L 



t 



no helicity , JJ e C 1 



(17) 



for t > Tq isc and where Tq 1sc ps rf (T^ /rd rag ) ~ L^a/En. Here, in contrast to the condition in the early Universe, 
a constant (in time) drag coefficient a has been assumed. Note that Eq. Ijl7|l indeed predicts slower magnetic field 
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FIG. 5: The evolution of the ratio of the kinetic and magnetic energy spectrum Ffe = _E fc m /_E™ ag for a maximal helical magnetic 
field in the turbulent regime. In this case equipartition (I\. « 1) is only established on very small scales. At the integral scale 
the kinetic energy is always much smaller than the magnetic energy. 



energy decay than it's counterpart Eq. i|10[l in the turbulent case, in particular, a longer relaxation time Tg lsc ^> tq 
and smaller decay slope 7 v i sc = 7turt>/2 for energy density. In Fig. [5] and Fig. ^] results of our numerical simulations 
of viscous non-helical MHD in the free-streaming regime are shown. For times longer than the relaxation time on the 
integral scale, small-scale power spectra are well described by power laws Ek cx £r of exponent (3 rs —2.0. This power- 
law is approximately consistent with a Reynolds number R el ~ const(Z) independent of scale I. Fig. ^] illustrates that 
Eq. I|17fl is a good approximation to the numerical simulations, though numerically simulated fields tend to decay 
somewhat slower than predicted, as observed in the sections before. 



B. Helical fields 



In the case of maximally helical fields one may use the constancy of helicity density in Eq. (|11|) to find 

dE E i 



dt H 2 a 



(18) 



yielding 



-1/3 , t vl/3 

E ~ E o[^f) L ~ L o[^) max. helicity , R e « 1 , (19) 



_V1SC / \ rr V1SC 

' o / V T o 



where Tq 1sc is as before. Results of our numerical simulations for this case may be found in Fig. 1101 and Fig. ^2 As 
in the the non-helical case these are consistent with a small-scale power law spectrum Ek oc k 13 with (3 sa —2.0. Note 
that, in contrast to before, agreement of Eq. (|19|l and the simulation seems excellent. 



IV. MHD WITH AMBIPOLAR DIFFUSION 



After recombination the Universe is only weakly ionized (i.e X e 1). Neutral particles, i.e. hydrogen atoms, do 
not respond to magnetic stresses and may therefore slip by the magnetic field lines, unless the scattering between 
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FIG. 6: Evolution of magnetic energy spectra in the turbulent regime for magnetic fields with initially maximal helicity. The 
spectral index of the energy spectra is n ~ 4. 



neutral and charged particles is rapid enough. To investigate if this is the case one has to consider the (here assumed 
incompressible) equations of MHD with a significant neutral component 

„ \ (V x B) x B 

+ Vi • V v; J = ^ ft a in (v< - v n ) (20) 

Qn (~^f + V " ' V V «^) = ~Sn «ni (v„ - Vj) (21) 

where g n , ft, v„, are matter density and velocity of neutrals and ions, respectively, and we will assume ft <C g n 
throughout. The momentum transfer rate due to neutral- ion collisions satisfy 

a ni = — a in w X e a in (22) 

Qn 

The equations of MHD are closed by including the induction equation 

<9B 

— = V x (v, x B) (23) 

for ions. The condition of tight coupling between ions and neutrals, i.e. = — v„ <^ may be derived from 
Eq. H21|) (noting that the first two terms are usually of the same magnitude) to be equivalent to 

^ Rj ^« Xe a in (24) 

One may show (cf. also [sH lifl Elj) self-consistently that in this limit the LHS of Eq. H20|l is negligible, leaving the 
ion- neutral drift velocity wb in the terminal velocity regime 

(V x B) x B 

v D = — : (25) 

47T ft a in 

Inserting this equality into Eq. I|21|l . and for v n v, where v is the center-of-mass velocity, one obtains the usual 
Euler equation Eq. QJ. The induction equation 123(1 is modified to include a dissipative term. Replacing Vj = vjj + v„ 
one finds 

;.Vx(vxB) + Vxff VxB » XB x B ) . (26) 
dt V 47T ft a in J 
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FIG. 7: The evolution of the magnetic energy in the turbulent regime for different initial energy spectra n, where Ek = 
k 3 |fofc| 2 oc k" . Here, the initial magnetic field is maximal helical. For comparison, also the theoretical damping law, E oc t -0 ' 67 , 
is shown. In contrast to non helical case, the damping law for a helical magnetic field is nearly independent of the spectral 
index n for n < 1. 



MHD of an ion-neutral mixture in the tightly coupled regime behaves thus as ordinary MHD with an additional dissi- 
pative term. The effect of this term may be estimated by defining an ambipolar Reynolds number as the comparison 
of the two terms on the RHS of Eq. I)26|l . i.e. 

vLa in ^ v L a in X e 

(VA 1 ) 1 V A 2 

where 

VA 1 w v A /Vx~e (28) 

is the Alfven propagation velocity in the ion-neutral weakly coupled limit, ft may be seen that the condition i? am b 1 
(assuming self-consistently v ~ va) automatically implies the tight-coupling condition Eq. I|25|l. It is thus evident 
that MHD with dissipation due to ambipolar diffusion in the tight-coupling regime may never become viscous due to 
this ambipolar "drag" [Hij]- In the language (JK098) appropriate to linear MHD this implies that overdamped modes 
proportional to the magnetic stresses do not exist. This is in stark contrast to MHD with fluid shear viscosity, or with 
momentum drag due to a homogeneous background component, where viscous MHD (i.e. overdamped modes) exist. 

These arguments assume the absence of other sources of dissipation. Consider, for example, shear viscosity due to 
neutral- neutral scattering as described by a term £> ra ^V 2 v„ on the RHS of Eq. I|21|) . Assuming viscous MHD due to 
this term (i.e. the kinetic Reynolds number Eq. J3J R e <§; 1) the condition for tight coupling is modified, and now 
reads rj/L 2 <C a ln X e . Nevertheless, even in this case one finds that the condition of tight coupling is equivalent to 
the requirement i? a mb 3> 1- The flow may thus be viscous in the tight coupling regime, but only due to sources of 
dissipation other than ambipolar diffusion. 

Once the flows reach the limit i? a mb ^ 1 the neutral species decouples from the flow. In this limit MHD evolution 
is described by Eq. 12U|) with v n — > and Eq. (|23l) completely analogous to MHD with free-streaming photons or 
neutrinos (cf. Sec. IHIjl and with a Reynolds number given by Eq. {TSJ. One may then show that the flow is viscous 
due to ambipolar drag. Only when a- ln is reduced by a further factor of ~ \/X e (or equivalently, the Alfven crossing 
time is reduced by the same factor) does turbulent MHD obtain again. When this happens typical fluid velocities are 
v » «a\ thus increased with respect to the tight-coupling regime. 

It is important to stress the following. The system of equations Eqs. (|2(jp. I|21|) and IL'-il) provides only a proper 
description of MHD in the fluid limit when particle species have mean free paths much smaller than the scale under 
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FIG. 8: Evolution of helicity as a function of time for different spectral indices n. In the numerically ideal case TC/TC„ 
should obtain. 
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FIG. 9: The evolution of the magnetic energy spectra in the viscous free streaming regime (R e <C 1) for a magnetic field without 
initial helicity. The simulations were performed on a mesh with 256 3 grid points, and the cut-off is k c ~ 16. 



consideration. Whereas for the scales we consider (assuming magnetic fields not too weak) this may be the case for 
protons and electrons due to Coulomb scattering, this condition gets violated at late times for neutral particles. In 
this limit, i.e. Z^ fp 3> L mixing of neutrals from different regions becomes significant. Higher moments of the particle 
distribution f n (x, v) (with the zeroth moment, density and the first, velocity) become significant such that a reduction 
of the Boltzmann equation to the Euler equation is not anymore adequate. When Z™ f 3> L and ZJ^ <C L one may 
nevertheless describe MHD by the fluid equations Eqs. (|20)l and iffi^) for ions, and to a good approximation, assume 
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FIG. 10: The evolution of the magnetic energy spectra in the viscous free streaming regime (R e <C 1) for a magnetic field with 
maximal helicity. The simulations were performed on a mesh with 256 3 grid points, and the cut-off is k c w 16. 

v n w due to mixing of neutrals from different regions. 

V. EVOLUTION OF COSMIC MAGNETIC FIELDS 

In this section we present detailed analytical results for the evolution of subhorizon magnetic fields between an epoch 
of magnetogenesis (e.g. the electroweak transition at T rts 100 GcV) and the much later onset of cosmic structure 
formation (at approximately redshift z fa 10). Our analysis draws on the general results found in the previous sections, 
but includes viscosities as applicable in the early Universe. In particular, we give coherence length and energy density 
as function of cosmic temperature, with generation epoch, magnetic spectral index, initial magnetic energy density, 
and helicity left as free parameters. Though the results arc fairly straightforward, when applied to the various regimes 
in the early Universe (i.e. turbulent and viscous due to photon- and neutrino- viscosity, respectively) a large number 
of expressions emerges. We therefore advise the more superficially interested reader to skip the third subsection of 
this section and proceed to the discussion of results in the next section. 

A. Initial Conditions 

We define the Fourier transform of the magnetic fields such that the spatial average of magnetic field strength may 
be written as follows 

< B{xf >= j A\nkB 2 k (29) 
The spectrum of the magnetic field at the magnetogenesis epoch is parametrized by 

B 9 c(lc) = B gc (L gc ) (J^j ' (30) 

Here, and in what follows, a subscript g denotes quantities at the magnetogenesis epoch and subscript c refers 
to comoving values. Here comoving lengths are defined as the lengths they would have at the present epoch (i.e. 
l c {T) = l(T) (a(T )/a(T)) where a is scale factor and T present cosmic temperature) and we define comoving field 
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FIG. 11: Evolution of the magnetic energy without (solid) and with maximal (dashed) initial helicity in the viscous free 
streaming regime (R e <C 1). The simulations were performed on a mesh with 128 3 grid points, the cut-off is k c ~ 16 and the 
spectral index n w 4. For comparison, also the theoretical expected damping laws are shown, i.e. iJmag oc t~ 6r (without 
helicity) and i? mag oc t~ 33 (with max. helicity). 



strength analogously as the field strength it would have at the present the epoch, if the field would only evolve 
according to the requirement of flux conservation (i.e. B C (T) = B(T) (a(T) / a(Tb)) 2 ) . Note that k = 2tt/1 and that 
whereas I denotes an arbitrary scale L always denotes the integral scale (k 1 the integral wave vector), i.e. the energy 
containing, scale. Our analysis will focus on blue spectra n > as appropriate for magnetic fields generated after an 
inflationary epoch by a causal process. Given these definitions one may compute the magnetic field energy density at 
an arbitrary epoch as a function of the temperature dependent integral scale and the scale factor 



, , B> c {L gc ) (31) 
a J mn y \ L gc J 

It is convenient to define a ratio r between magnetic energy density and a power of the total radiation entropy density 

r 55 w (32) 

since for constant comoving integral scale (i.e. no dynamic magnetic field evolution) this ratio stays constant with 
the expansion of the Universe. The dynamic (as opposed to geometric) evolution of the field is therefore more easily 
deduced from the evolution of r. The quantity r may be related to the ratio of magnetic field energy density and 
radiation energy density r e — qb/Qv by 

re = r 3(l5-J ^7 (33) 

with 35, g r denoting statistical weight in entropy and radiation, respectively. Note here, that r may be converted to 
average magnetic field strength 

B(T) = 5.72 x 10" 6 Gauss r^ 2 (T) 

' 3.909 ) \ 2.351 x 10- 4 evJ ^ 
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such that the comoving (present day) magnetic field strength is B c — 5.72 x 10 6 Gauss for r — 1. (For r 7 = qb/ = 1 
the comoving field strength of B c = 3.24 x 1CT 6 Gauss results.) The magnetic field strength given in Eq. (|34|l yields 
an Alfven velocity after the decoupling of photons of 

B , -cm ( r \ 1/2 [VLhh 2 Y 1/2 f T \ 1/2 

where gs = 3.909 has been assumed and with g — g^ in the fully ionized case before recombination as well as in the 
partially ionized case in the tightly coupled regime after recombination (cf. Sec. UVfl . Here gj, and f2t> denote baryonic 
density and fractional contribution to the critical density at the present epoch whereas h is the Hubble constant 
in units of 100 kms^ 1 Mpc~ 1 . We alert the reader to the distinction between decoupling of photons (i.e. Z 7 > L) 
(typically occurring before recombination) and recombination itself. Eq. (|35|l may be compared to the plasma speed 
of sound 



1/2 



where 7 = 10/3 and 2 for adiabatic and isothermal compression, respectively, and where we have neglected corrections 
due to the presence of helium. Note that below redshift z < 100 the baryon temperature falls more rapidly than the 
photon temperature, i.e. as Tb ~ a~ 2 . 

To determine the integral scale at the generation epoch, L g , we assume that turbulence pertains, such that L g is 
obtained by setting the Hubble rate at the generation epoch equal to the Alfven eddy turnover rate. This yields 



/ 2025 \ 1/6 M P1 _ 1/0 

L gc = L c {T g ) = y-^rj ivr g gsa ( ' 



1/3 

,V2 / T, 



L^O-WS^) 555^) (38) 



where Mpj = 1/VG ~ 1.22 x 10 19 GeV the Planck mass. In the above, the subscript denotes quantities evaluated 
at the present epoch. Finally, we parametrize initial helicity of the field by a dimensionless number h g 

H gc = h g H% ax with HL^^ C (L SC )L 9C (39) 



such that h g < 1. 



B. Evolution: The General Picture 



The evolution of a stochastic magnetic field in the early Universe is described by alternating epochs of turbulent 
MHD and viscous MHD. Here the latter epochs occur when viscosities due to neutrinos, or photons, become significant. 
Such a picture has already been established by JK098. With the assumed blue magnetic spectra, the gross features 
of magnetic field evolution are described by the growth of the integral scale. Following the arguments in the Sec. |H] 
the instantaneous integral scale is given by the equality between cosmic time and eddy turnover time at the scale L 

1 ~ V(L) ~ H(T) « -1 (40) 



icddy L p (T) tji 

holding equally for turbulent and viscous eras. In the above expression the subscript p denotes the proper (as opposed 
to comoving) value of the integral scale, v(L) is the the fluid velocity on scale L, and H is the Hubble constant. The 
velocities v may be determined from the Euler equation (Eq. l[T|l) by an approximate balance of either the second and 
third term on the LHS in the turbulent case (R e > 1) or the third term on the LHS and the dissipation term on the 
RHS in the viscous case. This yields 



in the turbulent case and 



v(L)^v A (L); R e >l (41) 



V (L)* V -^1±, v(L)* V -^l, R e <l (42) 
7? aL 



17 



in the photon (neutrino) diffusive and free-streaming viscous cases, respectively. Note that the velocities in the viscous 
case may also be written in a unifying way as v(L) f=a R e va with R e < 1 the Reynolds number Eq. evaluated with 
the Alfven velocity, i.e. v — va- 

Eq. (|4U[) is to be evaluated with proper quantities. Since a given scale expands continuously with the Universe, 
i.e. l p — al c , the eddy turnover (relaxation) time on this scale increases with the expansion of the Universe. This 
relaxation time increase is enhanced after the decoupling of photons by an additional decrease of the Alfven velocity, 
i.e. va oc B(T) I 'y/^KQh oc a -1 / 2 (whereas, va oc a when photons are still coupled to the MHD evolution). On the 
other hand, the Hubble time increases as tu oc a 2 during radiation domination (RD) and as oc a 3 / 2 during matter 
domination (MD). During turbulent evolution, the combined effect is such 

ioddy _ L/va a a/l/a 1 / 2 

w oc — oc l a (RD) oc -j- — oc a (MD) (43) 

tn tn a a ' 

that during RD larger and larger scales may be processed, i.e. that the comoving integral scale L c may grow as the 
Universe expands. In contrast, during MD the ratio between eddy- and Hubble- time stays constant, permitting only 
logarithmic growth of L c . This, however, is only the case while the fluid is turbulent. For sufficiently strong fields (see 
below) turbulence recommences right after recombination, with the fluid before recombination strongly dragged by 
free-streaming photons. In the viscous regime, with viscosity provided by photons, one finds ry 7 oc a 3 and a 7 oc a~ 4 
(cf. Appendix [UJ. This yields for the comparison of time scales 

cddy — — oc a (photon diffusion) oc a~ 5 / 2 (photon free-streaming) (44) 
tH tii 

where we have assumed radiation domination during the diffusive regime and matter domination during the free- 
streaming regime. It may be seen from Eq. (|44f) . that the integral scale may further increase during the viscous MHD 
regime when photons are free-streaming. On the other hand, one may show quite generally that an increase during 
the photon diffusion regime is always prohibited. Essentially identical conclusions result in the case of neutrinos. 

The following general picture for the evolution of the integral scale thus emerges. At early times, close to the epoch 
of magnetogenesis in the early Universe, the fluid is turbulent and as the Universe expands the comoving scale where 
one eddy turnover is possible in cosmic time is continuously increasing. By the process of a direct cascade the energy 
of this (integral) scale L c may thus be dissipated, leaving only the tail of the initial magnetic field at scales larger 
than the integral scale. The spectrum of the magnetic field is thus described by 

/ \ -n/2 



B e (!c) = B gc (L gc ) , l c > L c . (45) 

As the Universe cools down shear viscosity due to neutrinos becomes large, thereby reducing the Reynolds number 
of the flow. At the epoch when the Reynolds number becomes of order unity on the integral scale, a regime of viscous 
MHD commences. At this point, a further increase of L c is prohibited, since in the diffusive regime the relaxation time 
grows more rapid than the Hubble time. Any existing fluid flows are dissipated, leaving, nevertheless, the magnetic 
field at scales beyond i^OT hatact. Here i^OT re f ers to the integral scale when i? ej/ (L(TEOT)) has decreased to unity. 
This is analogous to the survival of magnetic fields in the overdamped regime of linearized modes, as discussed in 
JK098. Only some time after neutrinos have decoupled from the fluctuations, i.e. when l v ^> L, the integral scale 
may grow beyond that given at the epoch of end of turbulence, L^ OT . During this dissipation of magnetic fields 
in the viscous free-streaming regime, the integral scale grows more rapidly then during the turbulent regime. This 
more rapid increase (as opposed to a slower increase in the non-expanding case, cf. Sec. IIII|) is mainly due to the 
strong temperature dependence of the drag term. Since the neutrino drag is continuously decreasing, some time 
before neutrino decoupling at T w 2.6 MeV the fluid enters again a turbulent stage. At this point, the integral scale 
has grown to a value, as if the plasma would have not at all gone through a viscous period. The viscous period thus 
just delays the dissipation of magnetic fields. These evolutionary trends are shown in Figs. 1121 1131 1141 and 1151 which 
show the growth of magnetic coherence length, and the decay of magnetic energy density, for a number of initial 
conditions. The evolution of the kinetic Reynolds number R e is also shown for a particular scenario. The frozen-in 
state of magnetic fields during the diffusive neutrino regime with R e < 1 and the first part of the free-streaming 
neutrino regime becomes apparent by the plateaus in L c one finds at T ~ 10 7 — 10 8 eV. 

A similar picture results for magnetic field evolution after neutrino decoupling, but now with neutrinos replaced 
by photons. There are, however, subtle differences. After the electrons and positrons become non-relativistic, i.e. 
T < m e , the photon mean free path increases rapidly. This is particularly true during the period of the e^- annihilation, 
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i.e. 500 keV > T > 20keV j55|. Therefore, viscous MHD evolution with photons diffusing commences for a wide 
parameter space during this period. The epoch of viscous MHD with drag provided by free-streaming photons, which 
starts some time later, is always ended right at recombination (T i=s 0.26 eV). This, of course, is due to the virtually 
instantaneous decrease of a 7 by a factor ~ 10 4 which is the result of the loss of free electrons during recombination. In 
contrast to dissipation due to neutrinos, the viscous period due to photons thus does not only simply delay the growth 
of integral scale (i.e the dissipation). Rather, shortly before the end of the viscous MHD regime before recombination 
magnetic field strengths are genuinely larger and integral scales are smaller as compared to a scenario where the 
flow would have stayed turbulent all along. What happens after recombination depends then on the strength of the 
magnetic field at recombination. 

After recombination two different potential sources of dissipation come into play; ambipolar diffusion as well as shear 
viscosity due to hydrogen atoms. We assume, for the moment, a turbulent flow, i.e. v ~ va, with resulting conclusions 
turning out independent of this assumption. With the aid of Eq. I|40|l . evaluated shortly after recombination, one may 
determine the ambipolar Reynolds number Eq. (|27(l to be approximately i? a mb ~ 10 4 . It follows that the hydrogen 
atoms are tightly coupled to the flow (cf. Sec. II V|) . As this is the case, viscosity due to neutrals may play a role. 
One may evaluate the kinetic Reynolds number R e Eq. 0} on the integral scale at recombination due to hydrogen 
viscosity by noting that the Alfven crossing rate on the integral scale shortly before recombination (va/L) vcc is given by 
(^/a 7 H) rcc (cf. Eq. I|40l) and l|42l) applied in the viscous photon- freestreaming regime shortly before recombination). 
When this is done one finds that weak magnetic fields with small coherence lengths |5(l are entering a viscous regime 
due to hydrogen viscosity immediately after recombination, whereas strong fields do not. Here, the dividing magnetic 
field strength is given approximately by B% w 10~ 13 Gauss, corresponding to integral scales L c 10 pc (cf. Figs. 1121 
and llHJ) . Thus, fields with B c < B]! are not significantly processed immediately after recombination. Only some time 
later, when neutrals have decoupled from the flow, they are subject to further processing (i.e. increase of coherence 
length). This further increase in L mostly takes place at epochs with redshift z < 100. The increase in integral scale 
then occurs in a viscous regime with drag due to free-streaming hydrogen atoms (cf. Sec. IIV|I . quite analogous to 
the regime shortly before recombination. The flows become turbulent again only when the Universe is reionized as 
ambipolar drag then disappears |57j . The epoch of reionization occurs presumably at z ~ 10, at which point the 
integral scale grows virtually instantaneous to a larger value and stays approximately constant thereafter |58i |. 

In contrast, fields with strength B > B^ recommence turbulence after recombination. The growth of coherence 
scale during the recombination epoch is characterized by an almost instantaneous increase of a factor of order ~ 5, 
associated with dissipation into heat. Subsequent evolution does only increase the integral scale at best logarithmically, 
due to the peculiar red-shifting of Alfven crossing time and Hubble constant (cf. Eq. H43(0 . One may show that even 
such fields enter a viscous period later on, with viscosity first due to diffusing hydrogens in the tight coupling regime 
and later due to ion- hydrogen collisions (ambipolar drag) in the weak ion- hydrogen coupling limit. In any case, the 
magnetic coherence scale is not modified much anymore, even after the epoch of reionization [59j . Again, these trends 
may be followed in Figs. 1121-1151 

We have so far assumed that magnetic helicity is negligible. Due to the high Prandtl numbers in the early Universe 
helicity is conserved. Since for typical blue spectra with n > 1 non-helical fields decay more rapidly than helical fields 
(cf . Sec. [HJ , initial fields with sub-maximal helicity h g < 1 will ultimately reach a maximally helical configuration 
during the course of field dissipation. Somewhat oversimplifying it may be understood as the non-helical component 
dissipating leaving the fully helical component as a remnant. There exists a simple criterion of when the fully maximal 
case is reached. Using Eq. (|rt0l) and Eq. I|39(l and the assumption of initial non-helical evolution, one may show that 
maximal helicity is reached when the integral scale has grown to 

L™ = L gc h~^ . (46) 

The subsequent evolution of the field is different from the non-helical case. Whereas processing (i.e. cascading 
of energy to smaller scales) on the integral scale and the growth thereof still proceeds according to the requirement 
Eq. (|40|l . the required inverse cascade due to the conservation of helicity implies a transfer of energy from small scales 
to large scales. The instantaneous magnetic field spectrum is thus modified compared to Eq. I|30|l and given by 

B c {lc) = B gc (L gc )h 1 g / 2 (^^j 2 (J^j l c >L c (47) 

Note that, in contrast to Eq. Q45|l . the prefactor of this spectrum is time dependent (through the temperature 
dependence of L c ). In accordance with the findings of Sec.[H] the spectrum retains its initial slope n on scales l c > L c . 
It may be noted that due to the large dynamic increase of L C (T) between magnetogenesis and recombination, even 
fields with initially very small helicity typically have reached maximal helicity by recombination. 
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FIG. 12: The evolution of comoving coherence length for initial magnetic field configurations with different spectral indices n 
and inital magnetic helicities. Solid lines from top to bottom: (a) hg = 1, r g = 0.01, (b) hg = 10~ 3 , n = 3, r g = 0.01, (c) 
hg = 0, n = 3, r g = 0.01, (d) hg = 0, n = 3, r g = 10~ 5 . The labels Z„, Z 7 , Zjj refer to the comoving mean free paths of neutrinos 
and photons and the comoving Hubble length, respectively. The epoch of magnetogenesis was assumed to occur during the 
electroweak phase transition (T s = 100 GeV) . 
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FIG. 13: The evolution of the relative magnetic energy density r corresponding to the models shown in Fig. 1121 Solid lines 
from top to bottom: (a) hg = 1, r g = 0.01, (b) hg = 10 -3 , n — 3, r g = 0.01, (c) hg = 0, n = 3, r g = 0.01, (d) hg = 0, n = 3, 
r g = 10" 5 . The epoch of magnetogenesis was assumed to occur during the electroweak phase transition (T B = 100 GeV). 
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FIG. 14: The evolution of comoving coherence length for different initial magnetic field configurations. Solid lines from top to 
bottom: (a) h g = I, r s = 0.083, n = 3, (b) h B = 10~ 3 , r g = 0.083, n = 3, (c) h a = 0, r g = 0.083, n = 3. The labels i v , i 7 , « H 
refer to the comoving mean free paths of neutrinos and photons and the comoving Hubble length, respectively. The epoch of 
magnetogenesis was assumed to occur during the QCD phase transition (T s = 100 MeV). 
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FIG. 15: The evolution of the relative magnetic energy density r corresponding to the models shown in Fig. 1141 Solid lines 
from top to bottom: (a) h g = 1, r g = 0.083, n = 3, (b) h g — 10~ 3 , r g = 0.083, n — 3, (c) h g — 0, r g = 0.083, n — 3. The epoch 
of magnetogenesis was assumed to occur during the QCD phase transition (T g = 100 MeV). 
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C. Evolution: Analytic Results 

As discussed above, and observed in Figs. 1121 and 1141 the growth of the integral scale before recombination takes 
place in different regimes, i.e. turbulent MHD with neutrino viscosity, viscous MHD with free-streaming neutrinos, 
turbulent MHD with photon viscosity, and viscous MHD with free-streaming photons. Moreover, it is dependent on 
if maximal helicity has been reached or not. After recombination strong magnetic fields (B c > B'l) undergo only 
slight further evolution in the turbulent regime, whereas weaker magnetic fields may pass through an extended viscous 
hydrogen free-streaming regime. 

In the following we give analytic results for the integral scale and energy density in the different regimes, expressed 
as functions of the initial conditions. Here most (but not all) of the notation should be clear from the definitions in 
prior sections (e.g. subscripts of r, S, /, I, q, v, 7, b, and p indicate, total radiation, entropy, particles coupled to the 
fluid, leptons, quarks, neutrinos, photons, baryons, and protons, respectively, whereas a subscript denotes quantities 
at the present epoch). One may in principle also derive the transition temperatures, Teot and Teov, at which the 
fluid passes from a turbulent state into a viscous one, and vice versa, defined by R e « 1. We have nevertheless 
refrained from doing so, as the the number of expressions quickly exponentiates, and in some circumstances (i.e. 
0.5 MeV > T > 20 keV due to the annihilation) closed forms may not be derived. 



Evolution before Neutrino Decoupling 



The expressions for the integral scale and energy density during turbulent MHD evolution before neutrino decoupling 
are identical to those before recombination. The reader is thus referred to Eqs. (|55|l and (|56[) , for the case of 
submaximally helical fields, and to l|59|) and l|60[) . for maximally helical fields. The expressions for L c and r in the 
viscous neutrino free streaming regime, for submaximally helical (i.e. L c < L™ ax ) fields are 



Lc(T) 
r(T) 



L„ 



Gi 
Gi 



1 



G% M pl T3 
1 
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■ 2+7, 



with 
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49 tt' 



C(3) V 405 



5/2 



9}9u (gi+g g ) 



(48) 
(49) 

(50) 



Here L gc is given in Eq. I|38(l . As a numerical example, with index n = 3 and g r = 10.75, cjf = 5.5, g v = 5.25, gi = 3.5, 
and g q = as applicable between QCD-transition and neutrino decoupling one finds 



L C {T) ~ 2.1xl0- 2 pc (iiL) 
r(T) ~ 2.1 x 10~ 8 



1/2 



-3/5 



T 
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(51) 
(52) 



For maximally helical fields (i.e. L > L™ ax ) the appropriate expressions are found from Eq. (|48() and Eq. (|49ll by 
replacing n — > 1, except in Lg C which leaves a y/n dependence in L c , as well as multiplying the RHS of Eq. I|48(l by 



1/3 



_____ 2 /3 

and the RHS of Eq. I|49|) by h g ' . This yields the following numerical examples 
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(53) 
(54) 



where the same parameters as above have been assumed. 
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2. Evolution before Recombination 



The integral scale and magnetic energy during MHD turbulence R e > 1 (equally applicable before neutrino coupling 
and recombination) before maximal hclicity has been reached (i.e. for L c < L™ ax ) read 
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A numerical example for n = 3, g r = 3.36, gj = 2, and gs = 3.909 as applicable after the e ± -annihilation is given by 
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Similar, when maximal helicity has been reached L c > L™ ax (essentially using the above equations with n — ► 1 
replaced and inclusion of /i s factors as in the preceding section) one finds 



Lc(T) 



gs 



9° I 1/2 1/2 

\gr 9f , 
( 9s(T) \ 



1/2 1/2 
,9r 9f , 



2 / 3 , rr \ -2/3 

5-1 



h 2/3 



~ 2 / 3 / j, x 2/3 



(59) 
(60) 



A numerical example for g r — 3.36, = 2, and 55 = 3.909 is given by 



L C (T) ~ 4.4 x 10 _1 pc ( 
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(61) 
(62) 



which is virtually independent of the spectral index n. The expressions during viscous MHD with free streaming 
photons (where we assumed T < 20 keV as is usually the case) for L c < L™ ax are 



L C {T) — L gc 
r(T) = r. 
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and where ctt = 87r a 2 /3 to 2 w 6.65 x 10 25 cm 2 is the Thomson cross section. Here, R r = g r j (g r + qbm) accounts for a 
significant contribution of dark matter to the Hubble expansion shortly before recombination and X e is the ionization 
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FIG. 16: Shows the final magnetic field strength, B(To), for submaximal (i.e. L < £ max ) magnetic fields in the (n, T g ) 
parameter space. Here, the initial magnetic field strength is r g = 0.01. Results for different r g may be obtained by scaling the 
field strength by (r g /0.01) 1/2 . 



fraction (X e w 1 before recombination and X e w 10 4 after). As a numerical example for n = 3, fl^h 2 — 0.02, and 
X e = 1 one finds 

Lrj{T) K 4.0pc (i)" 2 ( I5 ^ v )" /S (^v)"' 5 (66) 

r(T) K 3 , , 10 -.» (i) (-j^y" 10 f^y*' 5 (_^y /5 (67) 

v ; VO.01/ V0-235y ViOOGeVy V0-259eVy v ' 

where T ~ 0.259 eV corresponds to the temperature at recombination and R r ~ 0.235 for f^ot^ 2 = 0.15. Similarly, 
for maximally helical fields L > i I 'j lax one finds the analytic expressions from those for the submaximal case, as above, 
by simply replacing in Eqs. (|63|l and I|tj4fl n — > 1, except in L gc , as well as adding a factor h^ 3 in Eq. I|tj3fl and a 

2/3 . j?m 

factor /i 9 in Eq. (|64|l . Numerical examples are given by 

it(T) K „. 3kpc ^ (i)" 3 (_|_) ,/s ( T5 a_)"' 3 (6 8 ) 

K ,2X10-" (Jfc) (i) 2 " ( ; J fe )" /6 f^)"' 3 (55^) (69) 



3. Evolution after Recombination 



In the turbulent regime after recombination the quantities of interest for L c < L™ ax are given by 

L ° = L - (vmtm 3S ° taW "'- ) ) 7 <S ) 7 (70) 

^7 ,../-./ /?o\^ 



r = r9l 7f^: 5soln(a/arcc V fe' (71) 
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FIG. 17: Shows the final coherence length, L(To), of submaximal magnetic fields (i.e. L < £ max ) ; n the (n, T g ) parameter 
space. Here, the initial magnetic field strength is r g = 0.01. Results for different r g may be obtained by scaling the coherence 
length by (r g /0.01) 1/2 



including a mild logarithmic growth factor ln(a/a rec ). Here J7 7 , f2b, and f2 t ot are the present day fractional contri- 
butions to the critical density of CMBR photons, baryons, and total matter, respectively. (In the derivation of this 
expression the contribution of radiation to the total density after recombination has been neglected. This induces 
about ~ 10% error in r and 5% in L c immediately at recombination for the values below, but is asymptotically 
correct). Numerical examples for n = 3, gso = 3.909, fitot^ 2 = 0.15, Sl b /i 2 = 0.02, fi 7 /i 2 = 2.48 x 10~ 5 are 

Lc(T) „ 12pc (J^f" (72, 

^.LlxlO--^^)"*' 5 (73) 

BAT) = 6.0 x 10-" Causa {^)~'" (74) 

where we have also evaluated the comoving field strength via Eq. (|34() . For the above quantities we neglected the 
factor ln(a/a rec ) (as for most fields the period of turbulent MHD after recombination is rather short) Similarly, for 
L > i,'j lax (when having attained maximal helicity) one finds 

L - = L " (vra •» '"(«/«"=») m ( I ) V3 h *' 3 < 75 > 

- g S0 Ha/a ICC )) ^ h 2 J 3 (76) 
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with numerical examples (with the input numerical values as above) given by 
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FIG. 18: Shows the minimal relative initial magnetic helicity, h g , necessary for a magnetic field to have become maximal helical 
(i.e. L > £ max ) a t the present epoch. For instance, initial magnetic fields generated at T g = 100 MeV with n = 3 become 
maximal helical if hg > 1.2 x 10 -7 . Note, this condition is independent of the initial magnetic field strength, r g . 

BJT) ~ 1.6 x 10- 11 Gauss ( Y'* ( ( Tg \ ^ (79) 
Note that there is only a residual dependence on spectral index n. 

VI. SUMMARY AND DISCUSSION 

The detailed numerical and analytical examination presented in the previous chapters has lead to a surprisingly 
simple picture concerning the gross features of cosmic magnetic field evolution, in particular, the evolution of magnetic 
coherence scale and energy density. The growth of the coherence scale is described by the simple causality relation 

v{L)/L^H{T) (80) 

independent if ocurring in high kinetic Reynolds number Re 1 turbulent flow with v ~ «a or during the multiple 
epochs of viscous (Re C 1) MHD evolution with w « % (with v given by Eq. I|42|l ) and independent of the helical 
properties of the fields. In particular, non-linear (direkt) cascading of magnetic energy to the dissipation scale always 
occurs on that scale. Remaining magnetic energy densities after evolution from very high temperature to an epoch 
with temperature T are then simply given, in the sub-maximal case, by all the initial magnetic energy present on scales 
I > L > , with L> the as yet largest length scale having been processed during prior evolution, and in the maximally 
helical case, by conservation of helicity density (see Eq. 1)11(11 with the field sitting on scales I > L > . Quantities of 
particular interest to cosmology are the anticipated present-day coherence length and magnetic field strength given 
particular initial conditions immediately after the epoch of magnetogenesis. We have shown that whereas strong 
magnetic fields B ~ 10 -13 Gauss are essentially undergoing no further evolution (i.e. growth of L > ) after 

recombination weak fields B < B^, do. In either case, after the later epoch of reionization the distinction between 
strong and weak disappears such that any reasonable strength field is turbulent at present. The present day field 
strength are then simply obtained by applying Eq. I|8U() with v Rj va today. This yields the correlation 

B ~5x 1CT 12 Gauss (J^ ( 81 ) 
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FIG. 19: Shows the final magnetic field strength, B(To), for magnetic fields which have become maximal helical (i.e. L > £ max ) 
in the (h g , T g ) parameter space. Here, the initial magnetic field strength is r g = 0.01, and for different r g results scale with 
(rg/0.01) 1 / 2 . Note the absence of a dependence on spectral index n. 








FIG. 20: Shows the final coherence length, L(To), for magnetic fields which became maximal helical (i.e. L > £ max ) j n t ne 
(h s , T g ) parameter space. Here, the initial magnetic field strength is r g = 0.01, and for different r g results may be rescaled by 
(rg/0.01) 1 / 2 . 
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FIG. 21: Exclusion plot on the comoving magnetic field strength on L — lOMpc (equally applicable to the present epoch and 
recombination) as a function of magnetic field spectral index. The solid lines shows the upper limit on B c (10Mpc) from an 
excessive distortion of the CMBR blackbody spectrum by magnetic field dissipation at redsfift z — 3 x 10 6 |2fi| . whereas the 
horizontal lines (from top to bottom) show the upper limits from present-day intergalactic Faraday rotation measurements 
and the possible overproduction of cluster magnetic fields (see text), respectively. 



between magnetic field strength and magnetic correlation length jf3(J . The magnetic correlation length itself is given 
by the initial conditions after the epoch of magnetogenesis (cf. Eq. (|7U[) and Eq. (|75|) for the sub-maximal helical- and 
maximally helical- case, respectively, with example values for an n = 3 spectrum given in Eq. I|72[l - Eq. (|74|1 and for 
fields which during the course of evolution have become maximally helical in Eq. (|77|l - Eq. (|79|l ). Both quantities 
are shown for non-helical fields as a function of spectral index and magnetogenesis temperature in Figs. lTHI and ll7l 

Concerning the evolution of maximally helical fields we have found and reverified prior work [lfij ] on the intriguing 
property of self-similar evolution of such fields. In particular, though maximal helical fields keep their initial spectral 
index on large scales (i.e. for L > L > ) the amplitude of the large-scale tail of magnetic fluctuations is subject to a 
seemingly "acausual" amplification on scales which are far larger than the distance an Alfven wave may travel (cf. 
Figs. Hfl and HOfl . The large-scale magnetic field spectrum for maximally helical fields is given by Eq. I|47fl . whereas 
Eq. (|45|l describes that of submaximally helical fields. We have numerically disputed the claimed effect (e.g. [19, 20] ) 
that maximally helical fields do not excite fluid motions and are therefore not subject to viscous damping (cf. Sec. lIIBjl . 
Rather, only due to the excitation of fluid motions the magnetic correlation length of maximally helical fields may 
continuously grow during the evolution of the early Universe. It is important to note that due to the large dynamic 
increase of L > between the epoch of magnetogenesis and the present (and the associated large dissipation of magnetic 
energy) a field with a minute amount of initial helicity typically evolves into a maximally helical field at present. This 
happens when the magnetic correlation length has grown beyond that given in Eq. (|46J) . Fig. 1181 shows the amount 
of initial helicity, h g , (cf. Eq. i|39|) ) as a function of magnetic spectral index required to reach a maximally helical 
state at present. Completely helical fields may thus not necessarily be considered an unlikely remnant of the early 
Universe. 

The evolution of magnetic fields during epochs with intermediate redshifts z 10 3 — 10 7 is described by turbulent 
evolution at higher redshifts, followed by a viscous MHD period without further growth of L y , and a viscous MHD 
period with comparatively rapid growth of L > with viscosity due to free-streaming photons (cf. Figs. IT2"1 and [14^1 . 
For essentially all interesting magnetic field strengths and spectra the plasma is in this last state shortly before 
recombination, allowing for the prediction of a correlation 



: 8 x KT 11 Gauss [ ^~ 
Vkpc, 



(82) 
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between magnetic field strength and correlation length shortly before recombination. Eq. (|82|l is obtained via Eq. 180|) 
and noting that v « i>a 2 / '(a 7 L) during viscous photon free-streaming. The correlation length to be employed in 
Eq. (|82|l may only be derived when the initial conditions shortly after the magnetogenesis scenario are known, i.e. via 
Eq. Hti3|) and the comments further below concerning maximally helical fields, whereas the instantaneous spectra are 
again given by Eqs. ijlS)) and (J37J. 

The correlation in Eq. (|82|l is almost identical to that one expects from a linear analysis (JK098). In contrast, for 
the formulation of limits on primordial magnetic fields due to magnetic field dissipation at redshift z w 2.5 x 10 6 and 
the concomitant production of spectral \x distortions in the CMBR, Ref. [2^| have employed the results of a linear 
analysis leading to the claim that fields of 3 x 10 -8 Gauss on scales of 400 pc are disallowed. Though the limiting field 
strength does not change when non-linear evolution is considered, as it is an energetic constraint, the comoving length 
scale does. This is due to the bulk of energy being not contained on the dissipation scale (at ~ 400 pc) but rather 
on the integral scale given by applying Eq. (|81|l at z « 2.5 x 10 6 . However, coincidentally the change is only mild, 
moving the limiting scale from 400 pc up to 1 kpc, since for such field strength the flow at z » 2.5 x 10 6 is only mildly 
turbulent (i.e. Re ~ 100) and since in the viscous regime both treatments almost coincide. When B > B]! the cosmic 
recombination process is associated with an almost instantaneous jump in the magnetic correlation length. How large 
this jump is than depends on the magnetic field spectral index determined during the magnetogenesis epoch. It would 
be interesting to examine at what field strength the associated energy dissipation could impact on the recombination 
process itself, especially in light of fields with strength B>6x 10 -11 Gauss being able to produce small-scale density 
perturbations as then v\ ^ v& (cf. Eqs. (jSHJ, and JSEl)- 

With the advances in high-precision CMBR anisotropy observations the interest in putative signals due to primordial 
magnetic fields has immensely risen. Essentially all current magnetic field induced CMBR anisotropy examinations 
assume fields of strengths ~ 10~ 10 — 10~ 9 Gauss on scales roughly the Silk scale, L w lOMpc. This is done, of course, 
since for much weaker fields the signal is hardly observable and when moving to much smaller scales not only are 
satellite missions like WMAP and Planck not able to resolve these but also are signals naively expected to be reduced 
due to the thickness of the last scattering surface a « 10 Mpc. By inspection of Eq. (|82|l it is clear that the scale 
of ~ 10 Mpc may not be the integral scale but rather a scale much beyond. We argue here that unless substantial 
primordial magnetic fields have their origin during an inflationary phase a search for primordial magnetic fields of 
~ 10~ 9 Gauss on 10 Mpc seems futile (at best controversial). This is due to a number of reasons. First, with fields 
which are causally generated during, for example, early cosmic phase transitions such strength on these scales are 
impossible to reach. This is of course due to the smallness of the Hubble scale in the very early Universe and since 
due to causality the spectrum must be sufficiently blue. Second, due to the blueness of the spectrum of causally 
generated magnetic fields, limits on smaller scales are easily violated. Fig. 1211 shows the maximum possible magnetic 
field strength on the scale 10 Mpc as a function of magnetic field spectral index. Fields which are above this value 
produce CMBR spectral [i distortions in excess of those observed |25j. It is seen that only for unrealistically small 
spectral indices B ~ 10 -9 Gauss on 10 Mpc may be reached. Finally, a direct constraint on the scale ~ 10 Mpc may 
be applied when the low-redshift collapse of a magnetized plasma to a cluster is considered 0] It is found that 
pre-cluster-collapse fields in of 4 x 10 -12 Gauss (corresponding to the authors 10~ 9 Gauss at redshift z = 15) are 
sufficient to reproduce observed Faraday rotation measures in present day clusters. Larger fields seem to overproduce 
the Faraday rotation measure and should therefore be ruled out. 

When trying to detect primordial magnetic fields a more promising and realistic alley should be the search for 
CMBR anisotropics on very small scales, in particular, on scales ~ 10 kpc (corresponding to multipoles I ~ 10 6 ), 
possibly close or only slightly above the integral scale, rather than the canonical 10 Mpc (multipoles I ~ 10 3 ). This 
is due to the existence of viable scenarios producing fields of interesting amplitude ~ 10 -9 Gauss on such scales and 
further such fields evading constraints from CMBR spectral distortions and observed cluster magnetic fields. Moreover, 
though this scale is much below the width of the last scattering surface, the expected signal is not necessarily small. 
In particular, for small scales I < 1 Mpc the magnetic field induced peculiar velocities are v ~ ua 2 / (a 7 0- CMBR 
temperature fluctuations follow ST/T cx v and for scales I < a are additionally suppressed by \flfo due to the 
thickness of the last scattering surface. Combing these factors one finds ST/T cx /-"- 1 / 2 for I > L, thus an increasing 
signal with decreasing scale. Here n is the spectral index of the primordial magnetic field. Primordial magnetic fields 
should therefore leave their strongest CMBR signal on small scales I ~ 10 6 . It remains to be seen if contamination of 
the primordial CMBR anisotropies by foregrounds may pose serious problems to such observations. 

Last but not least, we have already challenged in a prior publication [43l | the long-standing and wide-spread belief 
that cluster magnetic fields may not be entirely of primordial origin. It is typically argued that causal magnetogenesis 
scenarios (as, for example, due to local processes during the QCD- or electroweak transitions) yield only weak magnetic 
fields < 10~ 20 Gauss on the (pre-collapse) scale of a cluster of galaxies. Since during cluster collapse further magnetic 
field amplification by only modest factors 10 3 — 10 5 due to magnetic flux conservation result, starting from 10~ 20 Gauss 
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implies that one is still far from the observed /xGauss fields in clusters, requiring further field amplification processes 
such as a dynamo. The problem with this argument is that, a priori, it is not clear if the initial magnetic fields have 
to reside on the cluster scale itself, or if magnetic field energy density contained on much smaller scales may, during 
the collapse be transferred to the cluster scale. In fact, the only numerical simulation of the collapse of a magnetized 
plasma to a cluster to date |42| seems to indicate an approximate independence of the final result on initial magnetic 
coherence length, with final cluster Faraday rotation measures only dependent on the pre-collapse magnetic energy 
density (required are Bk4x 1CP 12 Gauss). It is currently not clear by what mechanism magnetic energy may inverse 
cascade from small scales to large scales during the cluster collapse. But if indeed it does, cluster magnetic fields could 
be entirely primordial, since magnetic fields of ~ 1CT 12 — 10 -11 Gauss on approximately ~ kpc scales are possible by 
either having magnetogenesis occur late, during the QCD phase transition, and/or magnetogenesis scenarios which 
generate a very small amount of initial helicity (cf. Eqs. (|7U|I . J7SJ), and l|81|l). It is interesting to note that such a 
scenario also lead to a "prediction" of magnetic field strength and amplitudes in voids, far from galaxies. Fields in 
such environments are presumably not affected by magnetic fields in galactic outflows and could be, in the optimistic 
case observable by future technology |3, 111 LU| • 
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APPENDIX A: THE MHD EQUATIONS IN MINKOWSKI SPACE 

We are using Gaussian natural units, i.e. fiQ = 1 = eo, and c = 1, for solving the MHD equations. On a static 
background in the Newtonian and non-relativistic limit the MHD equations are given by (see e.g. |44j 1 

|+V-M = (Al) 

<9v , 1„ Bx(VxB) 

— + v-V v = — Vp j '--V$ (A2) 

at g 47r g 

r)R 1 

^— = V x (v x B) + V 2 B (A3) 

at Ana 

A$ = Ang , (A4) 

where g, v, p, B, $, and a are matter density, fluid velocity, thermal pressure, magnetic field, gravitational potential, 
and electrical conductivity, respectively. The above equations have to be closed by an equation of state. 



APPENDIX B: CONFORMAL PROPERTIES OF THE MHD EQUATIONS IN AN EXPANDING 

UNIVERSE 

In the following we assemble the MHD equations in the FRW universe with the scale factor a and the Hubble 
parameter H for (a) relativistic MHD, i.e. when photons are still coupled to the plasma on the scale of the magnetic 
fluctuations, Z 7 <C L, and (b) for non-relativistic MHD in the opposite limit Ly 3> L. In both limits, the equation of 
state in the early Universe is well approximated by being isothermal due to incompressibility in the limit (a) and due 
to the efficiency of electron-photon Thomson scattering and the associated cooling in the limit (b). We then show 
how in both limits the MHD equations in the FRW background may be essentially reduced to those in Minkowski 
space, when appropriate scalings with scale factor of the physical quantities are introduced. For further details on the 
derivation of the equations we refer the reader to JK098. To lowest non-trivial order in 1/cr, the relativistic MHD 
equations are 

^ + I V . ((g + p) v) +3tf (q + p) = (Bl) 
(d 1 . _ \ v dp 1 Vp l/Bx(VxB)\ 

+ - ( v ' V) + i? V+ -£ + — + - — — ^ r 1 ) = 

\ot a J g + p ot a g + p a \ Air (g + p) J 
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^v 2 v+iv(V-v)) , (B2) 



\ ' ' ' 3 

f) \ 1 1 

— + 2H)B = -V x (vxB) + - ? V 2 B, (B3) 

dt J a iiraa^ 

where we assumed that g C m *C Q = £>fl u id (here g refers to internal energy density for radiation) and we kept only 
terms of the lowest order in v/c. The shear viscosity rj is given by |45| 

n= ^9tT 4 l mip /(g + p) , (B4) 

where g t is the number of relativistic degrees of freedom of the particles with the longest mean free path Z m f p . Using 
the following the rescaled variables (e.g. [Til Il5ll4^| ) : 

g = ga 4 p = pa 4 Be Ba 2 

v EE v T ee Ta fj ee r/a^ 1 (B5) 

dt ee dta^ 1 a ee aa 

the MHD equations (|B1|I - (|B3|I in the radiation dominated universe (i.e. a oc t" 1 / 2 and p — g/3) can be written as: 

^+V-((g + p)v)=0 (B6) 



d , \ v dp Vp BxfVxB 

~ + (v • V) 1 " 



dt J g + p dt g + P An {g + p) 

1 
3 



77 ( V 2 v+iv(V-v) I (B7) 



_=Vx(vxb) + — V 2 B, , (B8) 

A similar rescaling transformation can be done in the matter dominated (MD) regime (i.e. a cx i 3 / 2 and p<g) using 
super comoving variables |47l |: 

g ee ga 3 p ee pa 4 Bee Ba 2 

V EE Va 1 / 2 ?7 EE Tya- 1 / 2 dt EE dt G^ 2 (B9) 

ee a 3 / 2 H a ee aa 1 ^ 



The transformations yield almost the form of the ordinary non-relativistic MHD equations (cf. Eqs. ijAlfl - IA3|I 

Og 
dt 



d l+V-(gv) = , (BIO) 



dv 11 

— = + (v- V)v+-Vp+ B x (V x B) = -s , (Bll) 

dt g itrg 

dT\ 1 

^-Vx(vxB) = — V 2 B, . (B12) 
dt Ana 



Here, the dissipation term is 



-fj (V 2 v + (V • v) ] (B13) 



where g is again matter density. Here the term | H v in Eq. (|B13() represents the only difference to the MHD 
equations in Minkowski space and may be interpreted as a drag term. In particular, fluid momentum dissipation due 
to a homogeneous photon background with l-y >> L, i.e. photon drag, is described by the addition of a term —av on 
the RHS of Eq. i|B2l) due to free-streaming photons (and the dropping of the terms proportional to shear viscosity r] 
due to diffusing photons). In the scaled variables this leads to the following dissipation term 

-H +& ) v (B14) 



with a = a 3 / 2 a. 
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APPENDIX C: DISSIPATION OF ENERGY AND HELICITY 



With a homogeneous density and pressure distribution the total energy density is given by 

/•• /' 'I''', (V • v.,"') . (CM, 



v 



where va = B/ \J ^k(q + p) is the Alfven velocity. The magnetic helicity density is 



U = \- r [ d 3 xA-B, (C2) 
'v 



V 



where A is the vector potential, i.e. B = V x A. Using the MHD equations of Addx. IaI the time evolution of the 
above quantities is given by (up to surface terms) 

= -ib+P) X d3x(Vxv)2 "(4^)vX d3a;(VxB)2 ' (C3) 

"r-s//' 8 -^ 3 '- (C4) 

In the case when fluid momentum dissipation occurs by free-streaming particles rather than diffusing particles, i.e. 
X] V 2 v — > -«v, the first integrand (and prefactor) of the RHS of Eq. l|C3(l must be replaced by a gv 2 . Note that in 
the case of ideal MDH (i.e. a — > oo) the helicity 1|C2|) becomes a conserved quantity. 
In the radiation dominated regime the shear viscosity is given by (cf. also JK098) 

r\ = \ — Imfp, (C5) 

with Z m f p being the mean free path of neutrinos or photons with statistical weight gt, and where gf denotes the 
statistical weight of the total fluid energy radiation density. Here, the photon mean free path, Z^ fp = l/(<7Tn e ), as 
measured in comoving units, is 

^•""•^(tfff (rib)"*- <«» 

for T < 20keV. For the neutrino mean free path we assume l^ fp ~ l/(G 2 F T 2 (ni + n q )) (with the Fermi constant 
Gf ~ 1.1663 x 10~ 5 GeV~ 2 and the number densities n; and n q of relativistic leptons and quarks, respectively). We 
find at T — 2.6 MeV the comoving value 

'-^M^n^r- <c7) 

where gi = (7/8) 10 and g q = has been assumed. Note that when using Eq. i|C7ll for the computation of a v care 
has to be taken to not only include scattering but also neutrino annihilation (cf. |48j). On the other hand, neutrino 
self-scattering does not contribute such that, below the QCD-transition, gi should be effectively reduced to 3.5 (only 
e ± . We use the following drag coefficients for neutrinos and photons, respectively [jjl liif 

a, ^ 9 -j^-, (C8) 
9f W P 

a^\^'^\ Xe ^, (09) 

with ctt the Thomson cross section, m p the proton mass, and p 7 , Qb denoting photon- and baryon- density, respectively. 
In the high temperature regime (lMeV < T < mw) the electrical conductivity a is given by [H^l 

0.76T < cr < 6.7T , (CIO) 

where the larger value refers to the upper temperature bound. At temperatures below the electron mass the conduc- 
tivity becomes j44| 

an e T r m P /2 7 1 \ 3 / 2 

ff = -^^nfi :- . (OH) 



m e a In A V tt m 
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where r c is the mean time between two collisions, A = 67r i /2 ^172 (tp) (™~) ' anC ^ a ' TOfi ' anC ^ 716 ar ° s ^ ruc ^ ure 
constant, electron mass and electron density, respectively. The magnetic Prandtl number P m which gives the relative 
importance of the kinetic and magnetic diffusion is very large in the early universe, 



/keV\ 



3/2 



P m = 4tt 77 ct ~ 2.9 x 1 s ^— j (Cf 2) 

for T < m e . This allows to neglect the dissipation of magnetic field energy due to finite conductivity. We use the 
ion-neutral momentum transfer at low temperatures given by 

a in » 3.2 x f O- 9 s- 1 (j^) = am V h h 2 (£\ " (Cf 3) 

where a- m = 4.4 x 10 s (note a>i n ^ a n i), nn ~ rib denotes hydrogen density, and ao is the present day scale 
factor. The hydrogen mean free path l^ fp after recombination is determined by hydrogen-hydrogen elastic scattering. 
Scattering on electrons (or protons) may be neglected in computing l^ fp due to the small degree of ionization (i.e. 
X e «4x 10~ 4 ). Assuming a temperature independent |6fj |. cross section which is approximately ohh ~ t^b, wnere 
a-Q = 5.29 x 10~ 9 cm is the Bohr radius one finds 



~ M * w ' 3 »■= ( ira) " (s^v J < C14) 

for the comoving mean free path. Shear viscosity due to hydrogen-hydrogen elastic scattering may be estimated by 

^-I^-mxjo-s! (_^)-'(_|_y 5/2 (C15 , 



where vj 1 ? = y/3T/m p denotes hydrogen thermal velocity. Finally we estimate the mean free path of electrons in 
the plasma. After recombination, and with an effective cross section a ~ (a/T e ) 2 where a is the fine structure 
constant and T e the electron temperature we find 

^ fp , c ~i(r 2 pc (cie) 

between the epochs with redshift z ffOO and z fOO. Note that the comoving mean free path is independent 
of temperature. Below redshift z ss 100 the mean free path decreases even further due to a more rapid decrease in 
electron temperature than photon temperature. We take an constant ionization fraction after recombination of 

I e w4x f 0" 4 (Cf 7) 

neglecting residual dependencies on and 17. 

APPENDIX D: GENERATION OF HELICAL FIELDS 

To excite a stochastic magnetic field with or without initial helicity we choose a coordinate system in fc-space useful 
for helical fields with the orthogonal unit vectors {e + ,e_,k} (see e.g. 44]). By expanding the Fourier transformed 
vector potential A in this basis, i.e. 

A k = A+e+ + A k e + A^k (Dl) 
one obtains the magnetic field in the new basis 

B k = -ikx A k = -fc(A+e + - A k e_) . (D2) 
With this set of basis vectors the magnetic field spectra are the given by (cf. also [lfi| ) 

|B k | 2 = fc 2 (L4+| 2 + L4 k | 2 ) , (D3) 
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whereas the magnetic helicity becomes 

n = WrI d3fcA ^ k =(^F / d3kHk m 

with 

H k ^Ai-B k = -k(\A+\ 2 -\A k \ 2 ) . (D5) 

Note, that this choice of coordinate system reflects also that the helicity (|D4|I is a well defined physical quantity as it 
is gauge independent, i.e. independent of A^. To ensure a real vector potential A(x) and from that a real magnetic 
field B(x) the ^4^ have to fulfill the relation 

(At)* = -A± k ■ (D6) 
The magnetic helicity can be of either sign but the magnitude |i? k | is limited by the relation 

|flk| < k- 1 |B k | 2 . (D7) 

A magnetic field is said to be maximally helical if the equals sign in the above equation holds. From the relations 
(|D3(1 and (|D5(1 it can be seen that the strength of the magnetic field can be chosen independently of the magnetic 
helicity (in this approach one can consider cither (A^ ,A^) or (i? k , H^) as independent variables). This allows to 
excite stochastic magnetic fields with arbitrary helicity. To excite stochastic magnetic fields with a fractional helicity 
we choose 

K = \TfK > (D8) 

where / 6 [0,1]. This convention leads to 

|-ffk| = -j; |B k | 2 Y+J • (D9) 

for the magnitude of the helicity in terms of the magnetic field. From equation 1|D9(I it can be seen that (1 — /)/(l + /) 
is the fraction of the maximal helicity magnitude H max = fc _1 |B k | 2 . This can be used to adjust the magnetic helicity 
to an arbitrary magnitude. 

Note, that the choice of exciting magnetic fields with a fractional helicity (|D9|I is not unique. This particular choice 
just reduces the amplitude of the helicity spectra (compared to that of the maximal helicity spectra) by a factor of 
(1 — f)/(l + f). For non-maximal helicity it is also possible that the helicity spectra _ff k do not follow the spectra of the 
magnetic field |Bk| 2 , but are rather independently distributed in /c-space. The particular choice of the implementation 
of a fractional helicity may influence the evolution of magnetic fields. 



APPENDIX E: NUMERICAL METHODS 

We performed the numerical simulation using ZEUS3D j2|| I27I Ejf |6^|. All simulations were performed with 
periodic boundary conditions. This mimics an infinitely large volume, where the surface integrals of the MHD 
variables around the entire box vanish exactly. Furthermore, we extended the code for the purposes of our studies. 
We used Gaussian random field for the initial fluctuations of the magnetic components with zero mean. To ensure a 
divergence-free magnetic field this is done by exciting modes of the vector potential Ak in fc-space in the following 
way: The complex vector potential Ak = (A k , A 2 ., A k ) is generated by 

A\ = \A\\e^ ie [1,2,3] , (El) 
where the amplitudes \A l \ are randomly selected using a Gaussian distribution, i.e. 



^<7 k I 2 a 2 



PQAl\) = ^^exp{-L"\ , (E2) 
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grid: 256 3 
grid: 128 3 
grid: 64 3 
grid: 32 3 




0.1 
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FIG. 22: Time evolution of the magnetic helicity TL in the case of a maximal helical magnetic field. The loss of helicity is due 
to magnetic diffusion, which is solely due to numerical diffusion which can be seen by the resolution study. 



and the phases tp^ are randomly selected with an uniform distribution from the interval [0, 2ir]. The amplitudes are 
related to the variance <7k by 

|4| 2 oc of cx k n , (E3) 

where we assumed an isotropic universe, i.e. — Ai. These modes were excited up to a cut-off k c . The initial 
stochastic velocity field is generated in the same way as the initial magnetic field described above. In addition one can 
either generate the stochastic velocity field v by using equation l|E2|) directly, or one can generate a divergence- free 
velocity field by first exciting a vector potential A and then computing v = V x A. The latter avoids strong density 
perturbations. 
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